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Abstract. We consider the time evolution of a particle bound by an attractive one-dimensional delta- 
function potential (at x = 0) when a uniform electrostatic field (F) is applied. We explore explicit ex- 
pressions for the time-dependent wavefunction ippix^t) and the ionization probability V{t), respectively, 
in the weak-field limit. In doing so, ipriO^t) is a key element to their evaluation. We obtain a closed ex- 
pression for ipF{0,t) which is an excellent approximation of the exact result being a numerical solution 
of the Lippmann-Schwinger integral equation. The resulting probability density |V^F(0,t)|^, as a simple 
alternative to P(t), is also in good agreement to its counterpart from the exact one. In doing this, we also 
find a new and useful integral identity of the Airy function. 



PACS. 79.70.+q Field emission, ionization, evaporation, and desorption 
and ionization of atoms (including beam-foil excitation and ionization) 
single-electron tunneling 



34. 50. Fa Electronic excitation 
- 73.23.Hk Coulomb blockade; 



1 Introduction 

The ionization of atoms in an external electric field is one 
of the oldest issues in quantum mechanics. As a simple 
example, hydrogen-like atoms in a uniform electrostatic 
field have extensively been considered [Ij. Here, the back- 
ground potential caused by the field decreases without 
limit in one direction, and electrons initially in the bound 
state will eventually tunnel through the "barrier" created 
by the field, and will ionize. Accordingly, there are no true 
bound states. The tunneling rate for an ensemble of many 
independent electrons has been calculated based on the 
exponential decay law following from the experimentally 
supported statistical assumption that the tunneling rate 
is proportional to the number of available atoms. On the 
other hand, big experimental advances in the field of nano- 
scaled physics have increased importance of the study of 
the tunneling process of individual electrons; new nano- 
scaled devices have been devised, examples of which are 
tunnel junctions based on the electron resonant-tunneling 
eflFect [2] and molecular switches [3]. Therefore, a detailed 
understanding of the capabilities of these devices clearly 
requires a deeper knowledge of the tunneling process of a 
single electron subjected to an external field. However, in 
analytic studies of the time evolution of the tunneling pro- 
cess, e.g., leading to ionization, we have the mathematical 
diflftculty that there are no explicitly solvable models for 
a transition from a bound state to the continuum. Also, 
even obtaining the numerical solution with high accuracy 
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to this problem would not be an easy task either, especially 
in the strong-field limit where a highly oscillatory behav- 
ior is found in the time evolution of the bound-continuum 
transition. 

In this paper, we would like to study the time evolu- 
tion of a particle bound by an attractive one-dimensional 
delta-function potential (at x = 0) when an external field 
is applied. We will restrict our discussion below to a uni- 
form electrostatic field F for simplicity's sake. The single 
delta-function potential well (without an external field) 
has a single bound state. This would make easier an anal- 
ysis of the ionization. From an applied point of view, the 
delta-function potential system has actually been utilized 
heuristically to represent a short range atom or optically 
active defect as well as, in more complex combinations, 
resonant-tunneling junctions [4j^6j, molecular switches 
[6l[7l[8], and Dirac comb lattices [9]. This model of the field- 
induced time-dependent ionization was first discussed by 
Geltman [lOj. Later on, some various approaches to explic- 
itly obtaining the time-dependent solution ?/^i?(x,t) have 
been carried out |lllll2lll3lll4lll5lll6j . However, no exact 
solution in closed form (in terms of its actual calculability) 
has been found. The method we adopt here to solve the 
time-dependent Schrodinger equation is to turn it into an 
integral equation based on the Lippmann-Schwinger for- 
malism. The integral equation has thus far been focused 
mainly upon its numerical solvability. Also, we will employ 
both the X F interaction Hamiltonian (the scalar-potential 
gauge) and the p A (the vector-potential gauge) and then 
compare them (c/. for a detailed discussion of x F versus 
pA gauge problem, see Ref. fV7\). 
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A primary goal of this work is to pursue the analytical 
expressions for the wavefunction and the ioniza- 

tion probability V{t), respectively, in the weak-field limit 
{i.e., the field strength / ^ 1 relative to the strength of the 
potential well). As will be shown, i/jpiO^t) is a key ingre- 
dient to doing so. In [12] and [H], on the other hand, the 
ionization probabilities in the strong-field limit have been 
studied while neglecting the infiuence of the residual zero- 
range potential on the wavefunction ijjF{x^t) and using 
the numerical method in the Born-approximation scheme, 
respectively. They have been then compared to the exact 
result obtained from the numerical analysis. In addition, 
the ionization probability in the strong-field limit obtained 
in the scheme of the exponential decay approximation was 
shown to be a good approximation on the average to the 
exact result, although it cannot account for the short-time 
ripples found in the exact result [12]. In the weak-field 
limit, however, we cannot easily neglect the infiuence of 
the ultrathin potential well to the wavefunction. We will 
obtain a closed expression for ?/^i?(0,t) in this limit, and 
study its probability density [^^^(O,^)^ as a simple alter- 
native to V{t). 

The general layout of this paper is the following; in 
Sect. [2] we briefiy review the known results including the 
delta-function potential problem without an external field 
and the problem of a particle subjected to a uniform elec- 
trostatic field but not bound by the ultrathin potential 
well. Here, we also derive the Lippmann-Schwinger in- 
tegral equation for the current problem from the time- 
dependent Schrodinger equation, and review the time- 
dependent tunneling through the ultrathin potential well 
{or barrier) in the field-free case, which is an analytically 
solvable model [19j. In Sect. [31 analytical expressions for 
il)F{x^t) and V{t) are explored, and an explicit expression 
for V^f(0, t) in the weak-field limit is derived, partially with 
the aid of the exponential decay approximation. Section [H 
deals with |?/^i?(0, t)p, as a simple alternative to the ion- 
ization probability, which is in excellent agreement to the 
exact result obtained from the numerical analysis of the 
Lippmann-Schwinger equation; this contains the slack rip- 
ples in the time evolution, which cannot be found in its 
counterpart obtained entirely from the exponential decay 
approximation. Finally, in Sect. [5] we give the conclusion 
of this paper. 

2 General formulation 

The system under consideration is described by the Hamil- 
tonian 

^ = C - ^o^(:^) - xF{t), (1) 

where Vq > {) and F{t) = F ■ 0{t). We decompose this 

Hamiltonian into Hq = 2J1 — Vo ^{x) and Hf = |^ — f F, 

where /i = 2m and H = Hq -\- Hf- First, Hq has a single 
bound state, 

M^) = VBe-^\^\, (2) 

with eigen energy £^5 = — -^2^"' where B = Clearly, 
V;^(x) has a discontinuity at x = 0. All eigenstates and 



eigenvalues of ^0 , and also the completeness of the eigen- 
states, were discussed in detail in, e.g., Refs. [20j and [21j. 
The eigenfunction of Hf with (continuous) energy E is 
given by [22] 

= Ai{-(^)*(a:+f)}, (3) 

where the Airy function 

Ai(^) ■= ^ Coodye'^^ + '^y. (4) 

We have an alternative choice to x F in equation ([T]) for the 
field-interaction Hamiltonian; the Hamiltonian H in the 
scalar-potential gauge is then replaced by its counterpart 
in the vector-potential gauge, 

H, = - Vo S{x) , (5) 

where Hv,f = 27^ {P + Pc{t)}'^ with the vector poten- 
tial Pc{t) = Jq F{r)dr. The time-dependent Schrodinger 
equations in the scalar-potential and the vector-potential 
gauges read 

ih^iljF{x,t) = HiljF{x,t) ; 

ih ^ ipv{x, t) = Hy tlJv{x, t) , (6) 

respectively. Here, the two equations are connected by 
the relationship ^l)F{x^t) = e^^"^^*^^^ tljy{x^t). Remarkably 
enough, the ionization probability 7^ (r) = 1 — |(V;5|?/;f(t))P 
in the scalar-potential gauge is different from Vv{t) = 
1 — KV^b •e~'^^^*^'^^'^|?/;i?(r))p in the vector potential gauge. 
As was pointed out in [18j, for V{t) one switches on the 
field at t = and turns it off at t = r. Afterward the 
ionization probability is measured; for Vv{t) one turns 
off the vector potential instead of the field. However, it 
appears physically unrealistic to think of an experiment 
where the vector potential Pc{^) = Jq F{t') dr' is turned 
off. Therefore we will take V{t) under consideration. 

We intend to derive the Lippmann-Schwinger integral 
equation for this problem from the time-dependent Schrodinger 
equation. To do this, let us go ahead with the Schrodinger 
equation in the vector-potential gauge in the momentum 
representation, where the homogeneous solution (to the 
Schrodinger equation for Hv,f) can be obtained very eas- 
ily. By substituting 

^-(^' *) = 7^ f-00 dp V>v{p, t) ei-P (7) 

into the equation for tljy{x^t) in (|6]) and then multiplying 
_1_ e~i^P on both sides, followed by the integration over 
X, we acquire 

C^dp'V{p-p')Mp\t), (8) 

where 

V{p) = ^ C^dxe-iP-V{x), 

V{x) = JZ, dp eiP- V{p) = -Vo 5{x) . (9) 
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First, the (homogeneous) solution to the equation, Equivalently, we have the integral equation for the total 



in 9_ ^^^^(^^ ^) = t) easily appears as propagator 



^v,F{p,t) = Uy^rip^t) - (p{p,0) , (10) 



K:F{x,t\x',0) = KF{x,t\x',0) + X (18) 

where the time evolution factor , , , . 

U.Xp,t) = e-i/o^{^ + ^=(-)}^'^-,andtheinitial bound Here, the homogeneous solut^^^^^^ 

state <^F(a;,t) = J.^rfx'JiTFla;, tlx', 0)i/'6(a;') with 

Next, by inserting an ansatz (py{p^t) = (Pv,F{p^t) Xv{p^t) 

(19) 

with x«(p, 0) = 1 into equation ® , we have reduces to a closed expression [H] 

(12) +M(xe(i)-x;-ii?;At)} (20) 

Integrating equation m over t and then multiplying -^^ ^^^^^ Moshinsky function (Ml 

(Pv,F{p,t), we obtain 

(p„(p,t) = (p„,f(p,t) + i2^;7„,Hp,t) X V ; 2 \v2rtj y ' 

f dr {p t) f°° dp' ipy{p' t) (13) where erfc(2;) is the complementary error function. Also, 
iJ.F ' -oo ' j^^Yg equation (fTZl) 



/fF(a;,t|0,T) = Ko{x,t\0,T) ■eWH-^)-^it-^) } ^ 

(22) 

where the field-free propagator 



By using the Fourier transform in equation ([7|), we ar- 
rive at the the Lippmann-Schwinger equation for ^^^(x, t), 
which reads 

^l)v{x,t) = (t)v{x,t) + ^Vq /q dri^^,F(^,^|0,r) V^^(0,r)^.^ Ko(:^,^|0,r) = .^^^ ^ ^Ittt^^' ^ (23) 

Here, the homogeneous solution tx • i • ^ ^ -j xi 

' 1 PCX) , / N 1 -l-t IS also mterestmg to consider the momentum repre- 

^) = 72^ J-oo ^-,^(^' ^) e^"^ P^OP^- sentation of the propagator ]Cf{x, t\x\ 0) in equation 



gator i^^^F(^, ^|^^ t) appears as the Fourier transform of 

/7,,F(p,t)V-^(p,r) such that {p\U{W) = ]Cf{pAp' = (24) 



K^^F{x,t\x',T) = (15) 



^ ^r)-pi^^(^-^') .r/ -TJ-^ (n t\ = where /7(t)=e and = e ^^^/V^ . By apply- 

27r h J-oo ^ ^^,Fl/^'^; ing an iteration of replacing the integrand /CF(0,r|x^O) 

, g-^{5e(t)-5e(r)} . ^i2(T^ii^-^cit)}-{x-XciT)}f j^y ^hc cutirc expression of the right hand side and then 



27ri h (t-r) 



performing the Fourier transform to each term, we will ob- 

with the field induced translation ^^^^ explicit expression for ^i?(p, 0); after a fairly 

. . I rt . / N Ft^ ^ 1. 1^ • ^ A lengthy calculation we eventually arrive at the expression 

Xc\^) ~ Jq (^^Pc\^) ~ 2 ^■'^d tne nelQ-mQuceQ ac- 

tion Sc{t)'= dTpl{r{'= Then, it turns out JCF{p,t\p' ,0) = U^^f{p - Pc{t),t) Af{p,p' ,t) , (25) 

that0^(x,t) = i^^^i?(x, 0) ?/^(x', 0), where where 

'0(x', 0) = ipb{x'). On the right hand side of equation ([T4|) . ^ ^ 

this clearly represents free motion subjected to an exter- 

AFip^p'^t) = 5{p - pcit) - p') J^dnU:^^{p-pc{t),n)x 
nal field (so-called Volkov part), and the second term is r oo k-i 

the infiuence of the residual zero-range potential. J ^ ( 2:;;^^o)^ H lo' /-^oo ^^,f(Pz, r^) /7;^(pz, r^+i) 

Keeping in mind the gauge factor ei^'P''^^\ we then I k=i i=i 
easily find the propagator Kf{x^ t\x' ^r) in the x Uy f{p' ^'^k) - (^ 

scalar-potential gauge, from (p!5|) for the vector-potential 

gauge, as Here, the summation index k counts how many times the 

electron interacts with the delta-well while it travels from 
KF{x,t\x' ,t) = Ky F{x,t\x' ,t) ■ eH^^^c(t) -x'p^^r)} ^ p to p'. The Born term {k = 1) has a closed form 

and also the corresponding Lippmann-Schwinger equation "^f (p^P ' ^) ~ 2h^'^ \f^h f {p-p' -p^t)} ^ (^7) 

^l)F{x,t) = (t)F{x,t) + 1 1^0 Jo dr KF{x,t\{),T) ' ^l)F{^,r) . r / . \ / . M 

(17) {erfc (^) - erfc (v^t ^ %)} • 
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where a^W = 2i^{p-p' -Pc{t)} and bpit) =^^{(p- 
Pc{t))'^ — {p')'^}- This momentum representation /Cf(p, t\p\ 0) 
in the scalar-potential gauge in terms of {/7^,f(p, t)} from 
the vector-potential gauge rather straightforwardly demon- 
strates the complexity of the bound-continuum transition 
involving infinitely many complicated continuum-continuum 
transitions. 

Now, we briefiy consider the field-free {F = 0) case de- 
scribed by the Hamiltonian Hq. In this case, we can obtain 
an explicit expression for tljF=o{x^t) from the Lippmann- 
Schwinger equation (p!7|) [19]. We sketch this procedure; by 
applying the Laplace transform to this equation for x = 0, 
with the aid of the convolution theorem, we obtain 



where the dimensionless "time" r] = ^(^)^ {t — r) [26], 
and the dimensionless "energy" a. Therefore, we have, 
from equation 



• (32) 



KF{x,t\Q,T) = Ko{x,t\Q,T) e*^(*-^) 



d<y Ai(a) e 



^o(0,s) = 



(28) 



where the Laplace transformed iI)q{x^s) = C{il)Q{x^t)}^ 
and (^o(^7 s) = C {0o(^7 ^)}- From equation ([2Q|) for F = 
and X = 0, it follows that 



(/)o(0,t) = y^e-i^^' erfc^^-f^^t) , (29) 

which gives ^o(0,s) = — t ^ . s -- Applying the 

inverse Laplace transform to equation ([28|) [24j, we now 
acquire V^o(0, t) = e~^^^^, which subsequently yields, 
from equation (p!7|) . [25] 

(/)o(x,t) + {M(|x|;z5;^t) - M {\x\; -iB; ^t)} = 

VB {M {\x\;iB;^t) ^ M {-\x\',-iB',^t)} = 

M^) e-i^'K (30) 

Along the same line, even the total propagator /Co(x, t\x' ^ 0) 
for ^0 was exactly derived in Ref. [19j. On the other hand 
it is very non-trivial to apply the Laplace transform to 
equation ([T7|) for F 7^ [note the highly oscillatory factor 

e~^'^*~^^^ in equation ([22]) ] in order to derive the total 
propagator for the explicit expression of which has thus 
far not been known. In fact, equation (pT|) has been just 
numerically treated, namely, first it is numerically solved 
for '^^(O,^), which is substituted back into it for ipp^x^t). 



3 Explicit expressions for wavefunction and 
ionization probability 

In dealing with the the temporal integral in (pT|) with ([22]) 

analytically, the factor e~^^^~^'^^ in i^ir(x, t|0, r) is a 
major limiting one. We intend to circumvent this difficulty 
with the aid of the Fourier transform of Airy function 

Ai(<T) [2a, 



(31) 



By means of the eigenfunctions (/)£;(x = 0) in equation (|3|) 
and the substitution a = — { ^^^y }^ this becomes 

KF{x,t\0,T) = jr^dE Ko{x,t\0,T) ^ X 

gi(^-W^)(t--) 0^(0). (33) 

Here, we see that each a may be interpreted as a spectrum 
channel of the propagator KF{x,t\0,T). From equations 
(|T7)l and (|32j| it fohows that 



iPFix,t) = Mx,t) + J^^daAi{a)TFix,t,a), 

(34) 

where 

2 

TFix,t,a) := /pdr^^ e^^rl^ •ei^^(^''^)(*-^) •V'f(0,t) 

. (35) 

with the field-induced energy srixja) = ^ + f (^7^)^- 
We also note that the eigen energy E = — 2i£i?(0,cr). 
To explicitly evaluate the integration over a in ([34]) . we 
expand Ai(cr) in terms of the delta-function 6 {a) so that 



Ai(-) = E^(^r^(-)' 



(36) 



/c=0 



which follows from the definition of the Airy function in 
equation ([4|). By means of the relationship J da S'^^^ (a) h{a) = 
{-1)^ h^^\0), equation ([34|) becomes 



^ ' 27rm 



/c=0 



(37) 



We now calculate {ilJh\'^F{t)) for the ionization proba- 
bility V{t) = 1 - |(?/^6|V^F(t))p. From equations ([2|), ([35]) . 
and ([37]) . we obtain 



{M^F{t)) = {M<pF{t)) 



2 m 



Vb' X 



k=0 



cr=0 



(38) 



where 

GF{t, a) = /g* dT tpF{0, t) e*Ti(<^)-(*-^) x 

{e-T^(*-^)'erfc(-73 v^F^^ + 74 + 
-^)' erfc(73 Vt^^ + 74 VF^) } (39) 



of 2 (t- 
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with 7i (<j) 
and 74 



' 2m 



; 72 



BF 
2m 



; 73 



2 V 2ihm ' 



. Equation §9^) yields 



1=0 



T=0 



^3k+l+l (41) 



where 

grit, t) = 2 ^f{0, t) {M{-Xe{t - r); -iB; ^{t- r)) 
+ M{xe{t-Ty,-iB;^it-T))} (42) 

in terms of the Moshinsky function M{x; k;t). In obtain- 
ing equation (|4T]) from (|4Q|) . an iteration of the integra- 
tions by parts was applied. From equations (|38|) . (|4T]) and 
(|42|) we clearly see that the knowledge of i/jpiO^t) is the 
key to the actual evaluation of the ionization probability. 
Therefore, we restrict our discussion below to derive an 
explicit expression for ?/^F(0,t), a good approximation of 
which can actually be obtained in the weak-field limit. 

We will consider two different approximation schemes 
for '0i7(O,t) in the weak-field limit; first, we simply re- 
place KF{x^t\0^r) and ipFiO^r) on the right hand side of 
equation ([T7|) by Ko{x^t\O^T) and ?/;o(0,r), respectively 
[cf. equations ([23|) and (|3Q|) ]. which immediately yields for 
x = 



V^^(0,t) ^ (I^FiO^t) + VBe-i^^' erf (^-1^5 • 

(43) 

Clearly, the second term of the right hand side appeared 
from the interaction between the particle and the potential 
well in the limit of ^ 0. By replacing 0i?(O,t) here 
by ^o(0,t) as well and then approximately substituting 
this ?/^i?^o(0,r) into ?/^i?(0,r) in (|35|) . we will be able to 
consider a closed expression for in (|34|) . To this 

end, we mimic the procedure for the field-free result in 
(|3Q|) . Then, equation (|35|) is transformed into [25j 



2(3, 



£^rfc(f 



^«/^^/,^d.^erfc(f 



)} 



where s = \Jt — T\ ol{x) 



\x\ 



^J2ih/n 



(iFs) - 
and 



(3F{x,a) = ^ ^{\Ei,\ - eF{x,a)} . Now, from (H?]) and 
(|44l) with ^Pf{0, i - s^) ~ ^^(o, t - ^2)^ gj^^j ^-i^j^^. 



fe=0 



where the dimensionless quantity / : = 
the relative field strength, and 



cr=0 

(45) 
represents 



Here, (3f{x^ ^) ^ y y ~ ^^-^ ~ ^f^ • Clearly, equa- 
tion (|45|) with / = exactly reduces to its field-free coun- 
terpart in (|3Q|) . 

In the second approximation scheme, we make use of 
the exponential decay ansatz for V^f(0, r) in equation (p!7|) . 



(47) 



where the complex- valued energy E = Ef — ^ Ff with 
Ef = + ZA/; ZA/ is the level shift. Equation (|47j) with 
the substitution t = t — easily yields, from equations 
(fT7|) and ([22l) for x = 0, 



I 2ihB^t ^-iEt 
IT m 



e-i^'Yf{t), (48) 



tA/(0,t) = 0/(0, t) + 
or alternatively, 

V;^(0,t) = MM ■ (l - \/'-^yf{t)) ' , (49) 



where 17(t) = dz e'^' ' with ^ = 4^ 
^2 = It has been shown that in the weak-field limit 
I/I <C 1 the exponential decay law offers a good approx- 
imation of the exact result for the ionization probability 
T^{t) [TH fT4]: further, the semiclassical value Af^wKB = 
_ 5h B j2 -g excellent agreement to Z\/ up to / < 0.1, 
e 3/ a good approximation of // for 



8m 

and Ef^wKB 



/ ^ 1 [T2tll4y 27|. However, the exponential decay approx- 
imation clearly cannot account for any ripples in the time 
evolution of the ionization probability observed from its 
exact result. We would like to approximately recover the 
ripples by using equations (|48|) and (|49|) instead of (|47|) . 

After a lengthy calculation an explicit expression for 
Yf{t) reveals itself as 



Yfit) 



E 

k=0 



2fc + l 



1 - 



2! 5 



3 6 e 

55 



(50) 



(for the detailed derivation, see AppendixEJ, which allows 
us to obtain a closed expression for ?/^/(0,t) in equation 
(|48l) [or (|49])]. For = 0, equation ([50]) clearly reduces to 
the field-free result 



J^dz 



!(2/c + l) 



/c=0 



|^erf(v/C^). (51) 



Since there is no guarantee for each expression for V^/(0, t) 
in (|48|) and (|49|) to fulfill the normalization condition for 
the wavefunction, let us consider their combination for a 
correct numerical evaluation of ?/^/(0,t), namely. 



Tf{x,t,CT) 



L_{e-2«/3. erfc(^ -/?Hx,a) 
g2a/3, erfc(^ + l3f{x,(j)Vt)] . 



_^/(0,^)=c (^^/(0,t) + ^ 



2ihB^t ^-^Et 



(46) 



e-^^*F/(t)) +(52) 

1 



(1-c) 0/(0, i). (i- v/^^^/W 
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where the constant < c < 1 can be determined by nu- 
merical fitting to the exact result being a numerical solu- 
tion of the integral equation (pT|) for x = 0. As shown in 
figures. [Hand [21 equation ([52]) gives the excellent results, 
especially in the very short-time region when the field is 
switched on so that the decay process is slower than ex- 
ponential fT4]. 

Here, we also find a new and useful integral identity of 
the Airy function; with the aid of equations ([3T]) and ([5T]) 
we easily get 

Yfit) = C^daAiia) e" " (-^«')*) = (53) 
^ da Alia) erf (^6 - <t (-3^)^) i r • 

At the same time we have (5S] 

J\J i 12 '12 

= e-«^ {6^iiFi(l;f;6) + 1} , (54) 
where Whittaker's M-function 

M,,^(a) = e-* + l + 2/i;ei), and 

the confiuent hypergeometric function 

iFi (a; b-x) = ZZo ^ M- Equation p 

can easily be verified by means of the expansion e~^^ ^ = 
1 — -\- ■ ' ■ and then the comparison of both sides. 

Now, the comparison between ([53]) with ^2 = and ([54]) 
immediately yields the identity [30j 

Too^W erf(xV^) 

= 7^^-'' {faifi(i;f;a) +1}, (55) 

where = 

4 Time-Evolution of the probability density 
at the potential center 

It is still not easy to deal with a closed expression in ([38|) 
for the numerical evaluation of the ionization probability 
V{t) = l—\{ilJb\iljF{t))\'^- Considering the shape of ^^^(x) = 
y^g-^kl fQj. ^ large enough {or \f\ <C 1 for a given 
field strength F), it would be no harm to adopt for the 
qualitative study of V{t) the approximation 

lim V^e-^l^l S{x) = lim./^e-^^% (56) 

B^oo n^oo y TT 

which allows us to use 1 — |?/^/(0,t)p as a simple alter- 
native to V{t). See figure [3] for the numerical evaluation 
of |?/^/(0,t)p for various values of /; the closed expression 
in ([5^ with ([5Q|) gives results in good agreement to the 
exact ones [31J. Here we observe the slack ripples in the 
time evolution obtained from equation ([52|) . The origin of 
the ripples has been discussed in [12] and [32]; the initial 
bound state (pb{p) in equation ([TT]) is symmetric around 



p = 0. By applying the field, this symmetry breaks down 
in such a way that the motion of the particles in one di- 
rection is just accelerated so that they will easily leave 
the potential well. On the other hand, the other direc- 
tion gets slowed down until the particles stop, and then 
they reverse their direction of motion so that the particles 
again approach the potential well at x = 0. Then they are 
partially refiected from the potential well. This process is 
repeated until all particles will completely leave the po- 
tential well. From the results shown in figures [T][3] we may 
say that '0/(0, r) in ([52]) leads to a good approximation of 
the exact result for il)f{x^t) in the integral equation ([T7|) . 



5 Conclusion 



In summary, we studied the time evolution of a particle 
bound by an attractive one-dimensional delta-function po- 
tential (at X = 0) (or an ultrathin quantum well) when a 
uniform electrostatic field is applied. Thus far, no ana- 
lytically solvable model of field emission has been known. 
We obtained explicit expressions for the time-dependent 
wavefunction ipf{x^t) and the ionization (or bound state) 
probability V{t), respectively, in the weak-field limit, es- 
pecially that for V^/(0, t) [see equations ([48]) . ([49]) and their 
combination ([52]) with ([50]) ] which is a key element to the 
evaluation of ?/^/(x,t) and V{t). This explicit expression 
for ?/^/(0,t) was shown to be a much better approxima- 
tion of the exact result than its counterpart obtained from 
the exponential decay approximation [see equation ([47|) ] 
in that, e.g., the resulting probability density |?/^/(0,t)p 
as a simple alternative to V{t) can easily be numerically 
evaluated and is in excellent agreement to the exact result 
with the ripples, whereas no ripples can be observed from 
the exponential decay approximation. It is further sug- 
gested that even for the strong-field limit our result for 
?/;/(0,t) would be a better approximation on the average 
to the (highly oscillatory) exact one than the result from 
the exponential decay law which has been shown to be 
a good approximation on the average in the strong-field 
limit [12j. In studying this subject, we also found an in- 
teresting integral identity of the Airy function. Next, we 
will explore the analytical expressions for various (time- 
dependent) quantities in a delta-function potential system 
= Ea S{x - xx), where A = 0, 1, 2, • • • A/". 



Acknowledgments 



The author would like to thank G.J. lafrate for the stim- 
ulating discussion. 

* 



ILki Kim: An analytic study of the ionization from an ultrathin quantum well in a weak electrostatic field 



7 



A : A mathematical supplement - derivation 
of equation ( ISOl ) 

We would like to evaluate the integral dz e~^^ -^2z'^ ^ 
Ao(6,6) + ^1(6,6) + ^2(6,6), where 

00 

M^i,^2) = E WSl^ /o ■ (57) 

Along the same line with eq. ([54|) . we can obtain 

J^dze-^^^'z^- = e-«^i^%±fiil (58) 
with m = 1, 2, • • • , which immediately leads to 

^0(6,6) = e-«^ {1 + ^iFi (1; f ; 6) + 

V (-^^)'^ iFi(l;m+ I 
/ > (3m)! 6m + 1 ( ' 

m=l J 

Similarly, we get 

dz e-«^ + 2 = + , (60) 

and then 

CO / o \ 

= e 2^ (3n + i)! ^eTT+s ^- (^1) 

n=0 

Also, it follows that 
J^dze-^^^'z^ = 2e-«^{^ + ^iF,(l;f;a)} , 

Jo ^ ^ ~ 6m + 5 ' 

^2(6, 6) = e-^^ el {^^ + ^iFi (1; f ; 6) + 

/ V (3m + 2)! 6m + 5 f * v"^^' 

m=l J 

From equations ([59|) , (|6T]) , and (|62|) we easily arrive at the 
expression in ([5Q|) . 
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Figure captions 

Figure [H decay rates //(t) versus time where //(t) is 
obtained from the wavefunction ^^/(O,^) = e~^-^^'^^ with 
A(t) = {El, + Z\/(t) - ^Ffit)) t. Here the bound state en- 
ergy El) = — and we use h = m = B = 1 for plots. 



Two solid lines are si from the exact result and s2 from 
equation ([52|) while two dot lines are dl from equation 
(|43|) and d2 (straight line) from the exponential decay 
law. a) for / = 0.1, c = 1 for s2, and Ff = 0.0010 for d2 
; b) for / = 0.5, c = 0.65 for s2, and Ff = 0.1896 for d2 ; 
c) for / = 1, and c = 0.45 for s2, and Ff = 0.52916 for d2 
; d) for / = 2, c = 0.45 for s2, and Ff = 1.2115 for d2. 

Figure [2 level shifts Af{t) versus time t, where Af{t) 
is obtained from ijjf{0,t) as //(t) is. a) for / = 0.1 and 
Af = -0.0072 for d2 ; b) for / = 0.5 and zA/ = -0.0738 
for d2 ; c) for / = 1 and Z\/ = -0.10722 for d2 ; d) for 
/ = 2 and Af = —0.11235 for d2. Other parameters are 
the same as for figure [H 

Figure [31 |?/^/(0,t)|^ versus time t. All parameters are the 
same as for figures [T] and [21 Note that dl breaks down for 
/ = 0.5, 1, 2, i.e., \tjjf{0,t)f becomes greater than 1. 
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